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Abstract — In this contribution, the performance of a 
multi-user system is analyzed in the context of frequency 
selective fading channels. Using game theoretic tools, a 
useful framework is provided in order to determine the 
optimal power allocation when users know only their 
own channel (while perfect channel state information is 
assumed at the base station). We consider the reaUstic 
case of frequency selective channels for upUnk CDMA. 
This scenario illustrates the case of decentraUzed schemes, 
where Hmited information on the network is available at 
the terminal. Various receivers are considered, namely the 
Matched filter, the MMSE flher and the optimum flher. 
The goal of this paper is to derive simple expressions for 
the non-cooperative Nash equilibrium as the number of 
mobiles becomes large and the spreading length increases. 
To that end two asymptotic methodologies are combined. 
■ The first is asymptotic random matrix theory which allows 
I us to obtain expUcit expressions of the impact of all 
other mobiles on any given tagged mobile. The second 
is the theory of non-atomic games which computes good 
approximations of the Nash equiUbrium as the number of 
mobiles grows. Q 

I. Introduction 

Resource allocation is of major interest in the 
context of multi-user systems. In the uplink multi- 
user systems, it is important for users to transmit 
with enough power to achieve their requested qual- 
ity of service, but also to minimize the amount of 
interference caused to other users. Thus, an efficient 
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power allocation mechanism allows to prevent an 
excessive consumption of the limited ressources of 
the users. 

The most straightforward way to design a power 
allocation (PA) mechanism is as a centralized pro- 
cedure, with the base station receiving training 
sequences from the users and signaling back the op- 
timal power allocation for each user. Power control 
schemes in cellular systems were first introduced for 
TDMA/FDMA [1], [2]; more recently an optimal 
scheme was derived for Code Division Multiple 
Access (CDMA) [3]. In order to achieve the optimal 
capacity, the users may also be sorted according to 
some rule of precedence [4] . However, this involves 
a non negligible overhead and numerous non infor- 
mational transmissions. In addition, the complexity 
of centralized schemes increases drastically with the 
number of users. As discussed in [5], centralized 
algorithms generally do not have a practical use 
for real systems, but provide useful bounds on 
the performance that can be attained by distributed 
algorithms. 

A way to avoid the constraints of a centralized 
procedure is to implement a decentralized one where 
each user calculates its estimation of the optimal 
transmission power according to its local knowledge 
of the system. This is, for example, the case in 
ad-hoc networks applications. Most of the time, 
a distributed algorithm means an iterative version 
of a centralized one. Mobiles update their power 
allocation according to some rule based on the 
limited information they retrieve from the system. 
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Supposing that an optimal power allocation exists, 
a distributed iterative algorithm is derived from a 
differential equation in [6] and its convergence is 
proven analytically. A distributed version of the al- 
gorithm of [2] is presented in [7]. Building on these 
results, a general framework for power control in 
cellular systems is given in [8]. A review of different 
methods of centralized and distributed power control 
in CDMA systems is given in [5]. 

In this context, a natural framework is game 
theory, which studies competition (as well as co- 
operation) between independent actors. Tools of 
game theory have already been frequently used as 
a central framework for modeling competition and 
cooperation in networking, see for example [9] and 
references therein. Building on the framework of 
[8], a game theoretic approach was introduced in 
[10], [11]. Numerous works on power allocation 
games have followed since, a selection of which we 
present in Sec. UIl 

Game theory can be used to treat the case of 
any number of players. However, as the size of 
the system increases, the number of parameters 
increases drastically and it is difficult to gain insight 
on the expressions obtained. 

In order to obtain expressions depending only 
on few parameters, we consider the system in an 
asymptotic setting, letting both the number of users 
and the spreading factor tend to infinity with a fixed 
ratio. We use tools of random matrix theory [12] to 
analyze the system in this limit. Random matrix the- 
ory is a field of mathematical physics that has been 
recently applied to wireless communications to an- 
alyze various measures of interest such as capacity 
or Signal to Interference plus Noise Ratio (SINR). 
Interestingly, it enables to single out the main pa- 
rameters of interest that determine the performance 
in numerous models of communication systems with 
more or less involved models of attenuation [13], 
[14], [15], [16]. In addition, these asymptotic results 
provide good approximations for the practical finite 
size case, as shown by simulations. 

In the asymptotic regime, the non-cooperative 
game becomes a non-atomic one, in which the 
impact (through interference) of any single mobile 
on the performance of other mobiles is negligible. 
In the networking game context, the related solu- 
tion concept is often called Wardrop equilibrium 
[17]; it is often much easier to compute than the 
original Nash equilibrium [9], and yet, the former 



equilibrium is a good approximation for the latter, 
see details in [18]. In this paper, we derive the non- 
atomic equilibrium, which generally corresponds to 
a non- uniform for the users. 

The non-atomic Nash equilibrium is studied in 
this paper for several linear receivers, namely the 
matched filter and the MMSE filter, as well as non- 
linear filters, such as the successive interference 
cancellation (SIC) [19] version of those filters. 
However, in order to perform SIC, the users need 
to know their decoding order, in order to adjust 
their rates. In this paper, we introduce ways of 
obtaining an ordering of the users in a distributed 
manner. The ordering can be determined simply in 
a distributed manner under weak hypotheses. This 
gives rise to a different kind of power allocation, 
that depend explicitly on the order in which the 
users are decoded. 

Moreover, we quantify the gain of the non- 
uniform R\ with respect to uniform R\, according 
to the number of paths. The originality of the paper 
lies in the fact that we show that as the number 
of paths increases, the optimal R^ becomes more 
and more uniform due to the ergodic behavior of 
all the CDMA channels. This is reminiscent of 
an effect ("channel hardening") already revealed in 
MIMO [20] . The highest gain (in terms of utility) is 
obtained in the case of flat fading (which also favors 
dis-uniform power allocation between the users). 

The layout of this paper is the following. First, a 
detailed account of related works is made in Sec. |Ill 
In order to be self-contained, we introduce useful 
notations and concepts of random matrix theory in 
SecHni The communication model that will be used 
throughout the paper is detailed in Sec. UVl Asymp- 
totic SINR and capacity expressions are given in 
Sec. lYl The particular game played between users 
is introduced in Sec. |VIl along with the existence 
of a Nash equilibrium. Finally, theoretical results 
for the power allocation are derived in Sec. IVIII 
for unordered users and Sec. IVIIII when there is 
an ordering of the users. Analytical results are 
matched with simulations in Sec. |Kl Conclusions 
are provided in Sec. 

II. Related Work 

This section is dedicated to present some of the 
works that use game theory for power control. 
We remind that a Nash equilibrium is a stable 
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solution, where no player has an incentive to de- 
viate unilaterally, while a Pareto equilibrium is a 
cooperative dominating solution, where there is no 
way to improve the performance of a player without 
harming another one. Generally, both concepts do 
not coincide. Following the general presentation of 
power allocation games in [10], [11], an abundance 
of works can be found on the subject. 

In particular, the utility generally considered in 
those articles is justified in [21] where the author 
describes a widely applicable model "from first 
principles". Conditions under which the utility will 
allow to obtain non-trivial Nash equilibria (i.e., 
users actually transmit at the equilibrium) are de- 
rived. The utility consisting of throughput-to-power 
ratio (detailed in Sec. [Vll) is shown to satisfy these 
conditions. In addition, it possesses a propriety of 
reliability in the sense that the transmission occurs 
at non-negligible rates at the equilibrium. This kind 
of utility function had been introduced in previous 
works, with an economic leaning [22], [23]. 

Unfortunately, Nash equilibria often lead to in- 
efficient allocations, in the sense that higher rates 
(Pareto equilibria) could be obtained for all mo- 
biles if they cooperated. To alleviate this problem, 
in addition to the non-cooperative game setting, 
[23] introduces a pricing strategy to force users 
to transmit at a socially optimal rate. They obtain 
communication at Pareto equilibrium. 

In [24], defining the utility as advised in [21] 
as the ratio of the throughput to the transmission 
power, the authors obtain results of existence and 
unicity of a Nash equilibrium for a CDMA system. 
They extend this work to the case of multiple carri- 
ers in [25]. In particular, it is shown that users will 
select and only transmit over their best carrier. As 
far as the attenuation is concerned, the consideration 
is restricted to flat fading in [24] and in [25] (each 
carrier being flat fading in the latter). However, 
wireless transmissions generally suffer from the 
effect of multiple paths, thus becoming frequency- 
selective. The goal of this paper is to determine the 
influence of the number of paths (or the selectivity 
of the channel) on the performance of PA. 

This work is an extension of [24] in the case 
of frequency-selective fading, in the framework of 
multi-user systems. We do not consider multiple 
carriers, as in [25], and the results are very different 
to those obtained in that work. The extension is 
not trivial and involves advanced results on random 



matrices with non-equal variances due to Girko 
[26] whereas classical results rely on the work of 
Silverstein [27]. A part of this work was previously 
published as a conference paper [28]. 

Moreover, in addition to the linear filters studied 
in [24], we study the enhancements provided by the 
optimum and successive interference cancellation 
filters. 

III. Random Matrix Theory Notations and 
Concepts 

The following definitions and theorem can be 
found in [12] and will be used in the following 
sections. In this section, N and K are positive 
integers. 

Definition 1: Let z/ be a probability measure. The 
Stieltjes transform associated to v is given by 



Definition 2: Let v = [vi, . . . ,vn\ be a vector. 
Its empirical distribution is the function FJj : M 
[0, 1] defined by: 

FJ,{x) = ^4^{v,<x\i = l...N}. 

In other words, FJj{x) is the fraction of elements 
of V that are inferior or equal to x. In particular, if 
V is the vector of eigenvalues of a matrix V, FJj is 
called the empirical eigenvalue distribution of V. 

Definition 3: Let V he a N x K random matrix 
with independent columns and entries Vij. Denote by 
[■J the closest smaller integer. V is said to behave 
ergodically if, as N, K ^ oo with K/N a, for 
X G [0, 1], the empirical distribution of 



1 |2 


1 |2" 







converges almost surely to a non-random limit dis- 
tribution denoted -F^(-) and, for y E [0,a], the 
empirical distribution of 

\Vl,[yN\ I 5 • • • , \VN,[yN\ I 

converges almost surely to a non-random limit dis- 
tribution denoted F^(-). 

Definition 4: Let V he a N x K random matrix 
that behaves ergodically as in Def. [3l such as F^i-) 
and F^(-) have all their moments bounded. The 
two-dimensional channel profile of V is the function 
p^{x,y) : [0,1] X [0,a] R such that, if the 
random variable X is uniformly distributed in [0, 1], 
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then the distribution of p^(X, y) equals F^(-) and, 
if the random variable Y is uniformly distributed 
in [0,a], then the distribution of p^{x,Y) equals 

Theorem 1: Let Y = V0WbeaA^xi^ matrix, 
where is the Hadamard (element-wise) product 
and V and W are independent N x K random 
matrices. Assume that V behaves ergodically with 
channel profile p^{x,y) as in Def. |4] and that W 



has i.i.d. entries with zero mean and variance 4f. 



N- 



Then, as A^, — oo with K/N a, the empirical 
eigenvalue distribution of YY^ converges almost 
surely to a non-random limit distribution function 
whose Stieltjes transform is given by: 



m iz) 



lim —Trace ( (YY^ - zt) 



-1 



u{x, z)dx 

and z) satisfies the fixed point equation: 

1 



u{x, z) 



Jo 



p^{x,y)dy 



(1) 



p^(x',y)u{x',z)dx' 



— Z 



The solution to equation ([T]) exists and is unique 
in the class of functions u{x, z) > 0, analytic for 
Im(2;) > 0, and continuous on a; G [0, 1]. 

IV. Model 

We consider a single uplink multi-user sys- 
tem cell, i.e., inter-cell interference free case. The 
spreading length is denoted N. The number of users 
in the cell is K. The load is a = K/N. The general 
case of wide-band CDMA is considered where the 
signal transmitted by user k has complex envelope 

^kit) = ^ SknVkit - nT). 

n 

Vkit) is a weighted sum of elementary modulation 
pulses which satisfy the Nyquist criterion with re- 
spect to the chip interval (T = NTc): 

N 

Vk{t)=Y,'"ik'4'{t-{l-l)T,). 

1=1 

The signal is transmitted over a frequency selective 
channel with impulse response Cfc(r). Under the 
assumption of slowly-varying fading, the continuous 



time received signal y{t) at the base station has the 
form: 



K 

= 5Z XI ^fc" / 



Ck{T)vk{t -nT - T)dT + n{t) 



n k=l 



where n(t) is zero-mean complex white Gaussian 
noise with variance a^. The signal (after pulse 
matched filtering by t)) is sampled at the chip 
rate to get a discrete-time signal that has the form: 



K 



n 



(2) 



k=l 



where are A^ x A^ Toeplitz matrices representing 
the frequency selective fading for the A;-th user, is 
a A^ X 1 vector representing the spreading code of 
the /c-th user, and n is an A^ x 1 Additive White 
Gaussian Noise (AWGN) vector with covariance 
matrix (t^Iat. 

We consider the case of a multipath channel. 
Under the assumption that the number of paths from 
user k to the base station is given by L^, the model 
of the channel is given by 



Ck{r)= Y^r,km{T-Tk{^)). 



(3) 



where we assume that the channel is invariant during 
the time considered. In order to compare channels 
at the same signal to noise ratio, we constrain the 
distribution of the i.i.d. fading coefficients ?7fc(£) 
such as: 

E [r]km = and E = -f- (4) 

J^k 

Usually, fading coefficients rik{tj are supposed 
to be independent with decreasing variance as the 
delay increases. In all cases, q is the average 
power of the channel, such as E [|cfc(r)|^] = 
E [|r7fc(£)|^] = Q, for all channels consid- 
ered. For each user k, let hik be the Discrete 
Fourier Transform of the fading process Cfc(r). The 
frequency response of the channel at the receiver is 
given by: 

hk{f)= 5^r/,(£)e"^-2'^/^'=W|v[/(/)^ (5) 

1=0 

where we assume that the transmit filter '^{f) and 
the receive filter ^*(— /) are such that, given the 
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bandwidth W, 



if 



f </< 
otherwise. 



w 

2 



(6) 



Sampling at the various frequencies /i = — 



w I 

" 2 



N-1 



/2 - --f + jfW, - 2 ' TV 

obtain the coefficients hik, 1 < i < N, as 

Lk-l 



w 

2 ' 

we 



(7) 

Note that E [\hikf] = g. 

Since the users are supposed to be synchronized 
with the base station and for sake of simplicity, we 
will consider in all the following that users add a 
cyclic prefix of length equal to the channel impulse 
response length to their code sequenceH This case 
is similar to uplink MC-CDMA [30], [31]. As a 
consequence, matrices {0^} are circulant [32] and 
can all be diagonalized in the Fourier basis F [29]. 
Model dll) simplifies therefore to: 



PkSk + n 



(8) 



k=l 



where is a diagonal matrix with diagonal ele- 
ments {hik}i=i...N- For each user k, the coefficients 
hik are the discrete Fourier transform of the channel 
impulse response. 

We make the hypothesis that the users employ 
Gaussian i.i.d. codes with zero mean and variance 
[33]. This hypothesis enables us to state simply 
our results, however almost all of the results are 
valid for any distribution of the codes as long 
as it has mean zero and variance [16]. In 

particular, since every unitary tranformation of a 
Gaussian i.i.d. vector is a Gaussian i.i.d. vector (so 
that Wj = F^Vj has the same distribution as Vj for 
any i), we multiply y in ([8]) with F^ and obtain 
without any change in the statistics: 

K 

y = HfcWfeA/^Sfc + n 

k=l 

= (HyP0W)s + n (9) 
where is the Hadamard (element-wise) product. 

^Note that in the asymptotic case (when A'^ oo), the result holds 
without the need of a cyclic prefix as long as the channel is absolutely 
summable [29]. 



In (|9l), H is the frequency selective fading matrix, 
of size N X K: 



H 



hii hi2 



hiK 

hNK. 



^ is the root square of the diagonal power 
control matrix, of size K x K. 

W is an X i^' random spreading matrix: 



W = [wi|w2| • • ■ \wk] where 



Wik 



WNk. 



Note that asymptotically (as A^ oo), for a 
given multipath channel of length L, model dH) 
is also valid for the case of uplink DS-CDMA 
since all Toeplitz matrices can be asymptotically 
diagonalized in a Fourier Basis [29], [34]. 

In the following, we will assume that the fre- 
quency selective fading matrix H behaves er- 
godically, as in Def. [31 The two-dimensional 
channel profile of Hv^ is denoted p(/, x) = 
Pix)\h{f,x)\\ f e [0,1], X G [0,a]. f is the 
frequency index and x is the user index. This 
enables us to use Th.[T]in order to obtain expressions 
for the SINK. 

It is also assumed that the power of all users is 
upper bounded by Pmax and the square norm of the 
fading, on all paths, for all users, is upper bounded 

by /imax- 



V. Asymptotic SINR Expressions 



Let hi. be the A;-th column of H, and H 



i-k) 



be 



H with hfc removed. Similarly, let Wk be the A;-th 
column of W, and W(_fc) be W with removed. 
Let \/P(-fc) be with the k-th column and line 



removed. Finally, let G 



W 



(-k) 



H 



i-k) 



VPi-k) 



A. Matched Filter 

Supposing perfect CSI at the receiver, the 
matched filter for the fc-th user is given by = 
y/Pk{hk Wk). This leads to the following expres- 
sion for the SINR of user k 



SINR. 





gfgfc 


2 


^^gf gfc + gf 




-k)Gf_i^. 


gk 
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Proposition 1: [16] As A^, i^' oo with 

K/N — > a, the SINK of user k at the output of 
the matched fiker is given by 



SINK, = ^ 



where : [0,a] — M is given by 
/3'^^(x) = P{x)- 

m^)? 

am{x) + /; j', P{y) Hf, y)\' \ h{f, x)\' dfdy 

(10) 

and H{x) = J^'\h{f,x)\'df. 

Denoting SINR^ = Pf^, Prop. □ enables us to 
extract an approximation of the value of the SINR 
of user k in the finite size case 



MF 



-^fc AT Sn=l I ^nk I ^ 



N l^n=l\"'nk\ ~T~ 1^2 j^k l^n=l ^ i \""n.j\ I"" 



|2 , 1 



N 



nk I 



(11) 



We observe that Pi 



B. MMSE Filter 



Supposing perfect CSI at the receiver, 
the MMSE filter for the fc-th user is 
given by g^'^^^ = R"^gjt, where R = 

(^Hv^ w) (^Hv^ w) ^ + a^I^v^ • This 

leads to the following expression for the SINR of 
user A; [14] 

SINR, = gf (G(„fc)Gf_,) + a^I^)-' g,. (12) 

Proposition 2: [16] N,K oo with 

K/N a, the SINR of user k at the output of 
the MMSE receiver is given by: 



SINRfc = p- 



MMSE 



where p^^^^ : [0,a] 
the implicit equation 



is a function defined by 



P 



MMSE, 



XI 



P(x) 



m,x)fdf 



n n-2 -L P{ymf,yrdy 

~^ Jo l+/3MMSE(j^) 



. (13) 



Denoting SINR, = pf^^^, Prop. [2] enables us to 
extract an approximation of the value of the SINR 
of user k in the finite size case 



^MMSE 



1 ^ 

71=1 



" ^ N 2^jj^k 1 



Pj I ^71 j I 



_^MMSE 



From (|T2l) . we observe that 



dPl: 



(14) 



^MMSE 



From Prop. [2l we have the capacity of user k 



^MMSE 



The global capacity of the system is 



C 



MMSE 



riog2(l + /3^^^^(x))rfx. (15) 

^0 



C. Optimal Filter 

The term optimal filter designates a filter capa- 
ble of decoding the received signal at the bound 
given by Shannon's capacity. Hence it is difficult 
to define an SINR associated to it. However, results 
of random matrix theory can still be applied. Let 

Y = (b.VP W j . The definition of Shannon's 
capacity per dimension for our system is 



^OPT 



^ logo det (im + 4 YY^ 



N 



a. 



(16) 



As N,K ^ oo with K/N 

Cf^ ^ / log2 (l + i^idt) (17) 

where u is the empirical eigenvalue distribution 
of YY^, as in Def. [2l If we differentiate the 
asymptotic value C'^^ of (fTT]) with respect to a^, 
we obtain 



log2(e) 



-iy{dt) 



(1 + ^t) 

log2(e) ( / l^^^dt) j u{dt)^ 



log2(e) I m'^i-a') - - 



(18) 



where m'^(-) is the Stieltjes transform of the empir- 
ical eigenvalue distribution of YY^. From Th. [U 
m'^(-) is given by 



u{f, z)df 
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where u{f,z) is given by ([U) with p^^(f,x) = 

p{f^x) = P{x) \h{f,x)\\ Given that if 

fjopj ^ |g immediate to obtain from ^ 



as 



C" 



OPT 



log2 



m i—z 



—dz. 

z 



(19) 



Propositions: (7°^^ and C^^^e ^j.^ 
through the following equality 



C 



OPT 



C 



MMSE 



+ 



log2 



log2(e) 



1 + - 



MMSE/ 



MMSE( 



-dx \ df 



(20) 

Proo/.- See Appendix IXI-AI ■ 
The additional term in the right-hand side of (|20|) 
corresponds to the non-linear processing gain. It 
quantifies the gain in terms of capacity that can be 
achieved between pure linear MMSE and non-linear 
filtering. 

Assuming perfect cancellation of decoded users, 
successive interference cancellation with MMSE 
filter achieves the optimum capacity [35]. The fol- 
lowing proposition ensues from this fact. 

Proposition 4: [16] As A^, i^' 00 with 

K/N a, the optimal capacity is given by: 



OPT 



where jS^^*^ : [0,a] 
implicit equation 

(3^^C{x) = P{x) 



Tlog^ {l + /3''''ix))dx 
Jo 

M is a function defined by the 



\h{f,x)fdf 



• (21) 

Prop. |4J enables us to extract an expression that is 
analog to the SINR for the optimal filter. Similarly 
to the case of /^'^'^^^ in Sec. |V-B[ the derivative of 
this expression obeys the property -Pfc^- = /5f'^- 



VI. Games and Equilibria 

From now on, we denote SINR^ = whichever 
filter is actually used. 

A. Power Allocation Game 

A game with a unique strategy set for all users is 
defined by a triple {5, P, {uk)k^s\ where S is the 
set of players, P is the set of strategies, and {uk)kes 
is the set of utility functions, Uk : P''^' M. 



In our setting, the players are simply the users, 
indexed by the set S^^ = {1, . . . , K}. The strategy 
for a mobile is its power allocation P^, which we 
will assume belongs to a compact interval P = 
[0, Pmax] ^ IR- The utility measures the gain of a 
user as a result of the strategy this user plays. In 
[21], the author derives what he calls Throughput 
to Power Ratio (TPR) under minimal requirements. 
The utility of user k is expressed 



7k 



(22) 



We denote 7^ = 7(/5fc), where 7(-) is the same 
function for all users. In ((22l) . 7 is at least and 
should satisfy conditions detailed in [21] in order to 
obtain an "interesting" equilibrium. 

For example, in the simulations, we consider 
the goodput 7(/3fc), which is proportional to 
(1 — 6"^*=) where M is the number of bits trans- 
mitted in a CDMA packet. Remark that the usual 
definition of goodput would rather be considered 
proportional to q{(3k) = (1 — BER^t)^^, where BER 
is the bit error rate. However, this quantity is not 
zero when the transmitted power is zero. Using this 
function in the utility would lead to the unsatisfying 
conclusion that mobiles should not transmit at all, 
since the (improbable) event of a correct guess gives 
them infinite utility [10]. Therefore, an adapted 
version of the goodput is adopted, where a factor 2 
is added before the BER. The performance measure 
considered is hence proportional to q2iPk) = (1 — 
2BERfc)*^, leading to the expression above. This 
function has the desirable property ^2(0) = and its 
shape follows closely the shape of the original good- 
put g(-). This is a relevant performance measure, 
as each mobile wants to use its (limited) battery 
power to transmit the maximum possible amount of 
information. 

This utility is expressed in bits per joule. In the 
non-cooperative game setting, each user wants to 
selfishly maximize its utility. A Nash equilibrium is 
obtained when no user can benefit by unilaterally 
deviating from its strategy. 

To obtain the maximum utility achievable by user 
k, we differentiate with respect to the power P^ 
and equate to 0. We obtain 



'dP, 



■l\Pk)-l{Pk 



0. 



(23) 



For all filters under consideration, ([TOl) . ([131) and 
(|2T]) imply Pk^ = Pk, thus ([23]) reduces an 
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equation on /3k 



/3fc7'(A)-7(A 



0. 



(24) 



Eq. (|24l) is particularly interesting in the case 
when there exists a unique solution P*. 

The existence of a solution to (|24|) is guaranteed 
as long as the function 7(-) is a quasiconcave 
function of the SINR, i.e., there exists a point below 
which the function is non-decreasing, and above 
which the function is non-increasing [23], [21]. In 
addition, we assume that the function 7(-) takes 
value 7(0) = 0, so that users cannot achieve an 
infinite utility by not transmitting. This occurs for 
several functions 7(-) of interest, in particular the 
goodput [24], which we will use for simulations. 
Unfortunately, the capacity can not be used as a 
function 7(-), since it leads to the trivial result 
P* = for this utility function. The uniqueness of 
the solution (3* to is due to fact that the SINR 
of each user is a strictly increasing function of its 
transmit power. Given the target SINR P*, we obtain 
the strategy of users in the next section. 

VII. Power Allocation in the Nash 
Equilibrium 

A. Flat Fading 

In this subsection, we show that the results of 
[24] for Matched and MMSE filters are a special 
case of our setting when L = 1 (fiat fading case). 
In addition, we derive the power allocation for the 
Optimum filter. When there is only one path, for 
each user k, denoted by its index = x E [0, a], 
h{f,x) does not depend on /. Given the target SINR 
P*, we have explicit expressions of the power with 
which user k transmits for the various receivers. 

In Appendix IXI-B[ we show that the influence of 
the strategy of a player on the payoffs of other play- 
ers is (asymptotically) "small". It justifies the fact 
that we can obtain an equilibrium in the asymptotic 
setting, without the need for players to possess all 
the information on the system. Their local informa- 
tion is sufficient. In the asymptotic limit, we obtain 
results similar to Wardrop equilibrium: the strategy 
used by each user does not influence the strategy of 
other users. 

1) Matched filter: From Prop. [Tl the continuous 
formulation is 

P*{a' + J,-Piy)\hiy)\'dy) 



Pix) 



\h{x)[ 



or equivalently in a discrete form 

Pk = r^V^ -■ (25) 

Summing (|25l) over k = 1, . . . , K, we obtain a 
closed form expression for the minimum power with 
which user k transmits when using the matched filter 



1 



Pk = 9 for a < — . 

\hk\'l-aP^ P- 



(26) 



2) MMSE filter: From Prop. |2l the continuous 
formulation is 



Pix) 



\h{x)\' 

or equivalently in a discrete form 



Pk 



(27) 



Summing (|27l) over k = 1,...,K, we obtain a 
closed form expression for the minimum power with 
which user k transmits when using the MMSE filter 



Pk 



1 



a^P* 



a 



1+/3* 



for a < 1 + (28) 

P* 



Both (1261) and (128]) are the same results as in [24]. 

3) Optimum filter: Each user maximizes its util- 
ity for a SINR equal to P* . However, in the case of 
the optimum filter, the SINR is not defined directly. 
It is nevertheless possible to extract an equivalent 
quantity from the expression of the capacity, since 
the value of the capacity of user k at the equilibrium 
is given hy C* = log2 (1 + /5*). 

Proposition 5: The power allocation is given by 



Pk 



a^P' 



/3+ 



\hk\ l-aT+^ 
where P'^ is the solution to 



for a < 1 + — (29) 



a 



log2 + -alog2(e)- 



+log2 1 + 



1 



ap^ 



1 + 



a 



1+0^ 



alog2(l + /3' 



Proof: See Appendix IXI-C[ 



(30) 
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B. Frequency Selective Fading 

In the context of frequency selective fading, for 
each user k, denoted by its index jj^ = x & [0, a], 
there are L > 1 paths with respective attenua- 
tions hi(x), £ = 1,...,L, which are i.i.d. ran- 
dom variables with some known distribution. We 
suppose that has mean zero, and the dis- 

tributions of the real part and imaginary part of 
hi(x) are even functions, as for example the Gaus- 
sian distribution, which we consider in the simu- 
lations. h(f, x) depends on / through h{f, x) = 
ELi^Ka;)e-2"^^(^-i). Given the target SINK f3\ 
the Nash equilibrium power allocation is determined 
by implicit equations for the various receivers. 

1) Matched filter: The continuous formulation is 



It follows that asymptotically n 



(and 



simulations in Sec.|K]prove that this approximation 
is valid for moderate finite values of N). From (|3T1) . 
we obtain a formula similar to (|26|) 

1^ for „ < 1 . (34) 



Pk 



2) MMSE filter: The continuous formulation is 



Fix) 



/3* 



\h{f,x)\'df 



or equivalently in a discrete form 



Pk 



"T 1 



(35) 



(36) 



Fix) 



a^g(x) + /; /; F{y) \h{f, y)\' \h{f, x)\' dfdy 
{H{x)f 

or equivalently in a discrete form 



In m, Kk = h{^,-^). 

As previously, when the number of users tends 
to infinity, J2f=i Pj l^njf' is asymptotically a 
constant (not depending on n), denoted fi. 



Pk = r ■ 

N Z^n=l \'^nk\ "T ^ Z^n=l \"'nk\ ^ Z^j^k^j I"" 



1 - 



a/3* 1 |fenj| 
1+13* K 



N 

n=l 



'"nk 



(31) 



In (ED, /i„fc = /i _(^,#). 

In this expression, the power allocation of user k 
seems to depend on the power allocation and fading 
realization of all the other users. However, when the 
number of users tends to infinity, the strategy of any 
single user does not have any influence on the payoff 
of user k, as shown in Appendix IXI-B[ Hence, 
the appropriate framework is non-atomic games. 
The expression j^Yl!j=iPj \hnj\^ is asymptotically 
a constant (not depending on n), denoted VL. 



where E, = ^ELJ/^™.|. 
It follows that asymptotically Vt 

obtain a formula similar to (|28l) 

1 a^P* 





1-Q 



1+/3* 



Ekl- 



a 



13* 



1+/3* 



for a < 1 + — . 

13" 



(37) 



we 



(38) 



1 - Ef=i %^ 



(32) 



3) Optimum filter: Each user maximizes its util- 
ity for a SINR equal to (3*. However, in the case of 
the optimum filter, the SINR is not defined directly. 
It is nevertheless possible to extract an equivalent 
quantity from the expression of the capacity, since 
the value of the capacity of user k at the equilibrium 
is given by C'' = ^ log2 (1 + (3*). 

Proposition 6: Asymptotically, as N,K — ^ oo, 
the power allocation is given by 



Pk 



where 



^Af I, |2 



a 



As — * oo, we can apply the Central Limit 
Theorem to the sum of random variables 



/3+ 
l+(3^ 



for a < 1 + — (39) 



1 ^ 



\h 



nj I 



E, 



(33) 



where (3^ is the solution to 

a log2 (l + - a logs (e) ^ + ^+ 
1 a/3+ 



+log2 1 + 



It tends to its expectation, which is equal to 1 (see 
Appendix IXI-DD . 



/3+ 



a log2(l + /?*). 

(40) 
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Proof: The proof is similar to the proof of 

Prop, m ■ 

We observe that for all filters considered, the 
optimal PA is a constant times the inverse of the 
total energy of the channel Ej. Via Parseval's The- 
orem, Ej = J2e=i l^e (if) I • It is a sum of i.i.d. 
random variables. As the number of paths increases, 
the optimal PA tends to a uniform PA. This is 
an effect similar to "channel hardening" [20]: as 
the number of paths increases, the variance of the 
distribution of the channel energy decreases and 
the Nash equilibrium PA becomes more and more 
uniform for all users. 

VIII. Successive Interference 
Cancellation 

The optimal filter gives a bound on the perfor- 
mance that can be achieved through (non-linear) 
filtering at the base station. In order to improve the 
performance of the system, we introduce Successive 
Interference Cancellation (SIC) [19] at the base 
station. Under the assumption of perfect decoding, 
SIC improves immensely the performance of lin- 
ear filters (Matched Filter or MMSE Filter). The 
MMSE SIC filter actually achieves the optimum fil- 
ter bound, under the assumption of perfect decoding. 
The principle of SIC receivers is quite simple: users 
are ordered and are decoded successively. At each 
step, supposing that the user has been encoded at the 
appropriate decoding rate, the signal is decoded and 
its contribution to the interference is then perfectly 
subtracted. This removes some of the inter-user 
interference and therefore increases the SINR of the 
following decoded users. 

The challenge is that the users must transmit 
at the appropriate rate to avoid the catastrophic 
occurrence of imperfect decoding. Usually, the or- 
dering of users is done in a centralized way, at the 
base station which then advertises it to the users. 
However, for the protocol to remain distributed, 
users should be able to decide, based on their local 
information, at which rate to transmit. 

At equilibrium, the rate is determined by the 
SINR and it is the transmission power of the 
user that is determined according to its rank of 
decoding. The equilibrium PA can be determined 
in a simple manner when the number of multipaths 
is finite (L < oo) and the number of users is very 
high (K oo). In Sec. IVIII-A[ we make use of 



the fact that the whole law of Ej is realized in this 
case, so that users automatically know their rank 
of decoding. Another manner to give a (random) 
ordering of decoding is to introduce an additional 
degree of liberty in the system. In Sec. IVIII-B[ we 
develop a correlated game framework that enables 
users to learn their rank of decoding in a simple 
way. In the following, we assume that each user 
has a unique has a unique i.d. number j ranging 
between I to K. 

A. Ordering when K ^ oo 

If the number of users K ^ oo, with L fixed, 
the whole law of the total channel energy will be 
realized. Assume the base station advertises to the 
users that they will be decoded by decreasing total 
channel energy. Each user knows, according to the 
realization of its fading, its rank in the decoding 
order given by K times 1 minus the cumulative dis- 
tribution function D{-) of the total channel energy 
E,. 

rankj = K{1- D{Ej)). 

In case that the base station advertises to the users 
that they will be decoded by increasing total channel 
energy, user j will have rank rankj = KD(Ej). 

B. Correlated Equilibrium 

We wish to introduce a simple mechanism that 
enables players to coordinate and to know in which 
order they will be decoded. We place ourselves 
in the context of correlated games. The notion of 
correlated equilibrium was introduced by R. Au- 
manrHin [36] and further studied in [37], [38], [39]. 
They represent a generalization of Nash equilibrium. 
The important feature of correlated games is the 
presence of an arbitrator. An arbitrator needs not 
have any intelligence or knowledge of the game, 
it needs only to send random (private or public) 
signals to the players that are independent of all 
other data in the game. In the context of non- 
cooperative games, each player has the possibility 
not to consider the signal(s) it receives. Coordina- 
tion between players turns out to be useful also in 
the case of cooperative optimization. The signals 
enable joint randomization between the strategies 
of the players, possibly resulting in equilibria with 

'Prof. R. Aumann has received in 2005 the Nobel prize in economy 
for his contributions to game theory, together with Thomas Schelling. 
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higher payoffs. The concept of correlated games 
was recently introduced in a networking context in 
[40], where the authors consider a simple ALOHA 
setting. 

The simplest and most intuitive coordination 
mechanism is given by a common signal which 
users as well as the base station overhear before 
each transmission. There are K\ possible permu- 
tations of K users. Hence, the arbitrator broad- 
casts a signal to the users belonging to the set 
{0, . . . , i^'! — 1}. Each of these numbers corresponds 
to a permutation tt of {!,..., K} that gives the 
(random) ordering of decoding as rankj = 7r(j). The 
users can then adjust their transmit power according 
to this ordering. In terms of size of the message, 
this is equivalent to the case when the base station 
decides the decoding order and broadcasts it to the 
users, or sends K individual messages of In(i^) 
bits containing the rank, since ln(K\) = KId.{K) + 
o(K \n(K)). However, there is no need of either any 
knowledge of the system or computations at the base 
station in the case of the correlated mechanism. 



C. SIC Power Allocations 

In both cases, once the users know their order, 
they can calculate their transmit power according to 
the filter that is used. The equilibrium still occurs 
when all users reach the SINR P*. A single user will 
not benefit by deviating, since it would decrease its 
utility. From now on, index k denotes the rank of 
decoding. 

In the case of the matched filter with SIC, the 
SINR of the user decoded at rank k is 



(41) 

From (|4T]) . we get the equilibrium FA of user k as 

Pk = 

N Z^n=l \"'nk\ "T ^2 j>k /-in=l ^ 3 |"-nfc| 



( L^r^ \h I 

I TV Z^n=l \"'nk\ 



(42) 



In the case of the MMSE filter with SIC, the SINR 
of the user decoded at rank k is 



MMSE 



1 ^ 



1 



n=l 



2 _i_ J_ -Pjl^^iP 



a 



(43) 

From (l43l) . we get the equilibrium of user k as 



(44) 



For flat fading, a simple recursion gives the 
equilibrium R^ (see Appendix IXI-E|) . We obtain 
respectively 



pMMSE 
^k 



\hk 



K-k 



(46) 



As far as frequency-selective fading is concerned, 
this gives us the form of the asymptotic expressions. 
Asymptotically, the power allocation of one user 
will not depend on the R\ of the other users, as 
shown in Appendix IXI-B[ With a similar reasoning 
as in Sec. IVIIl the expressions mimic (|45l) and (|46|) 
with the total channel energy Ek replacing \hk\ , 
i.e., 

^2 p. / . ^K-k 



pMF 
^k 



Eh. 



El. 



pMMSE 
^k 



1 + 



1 



K-k 



(47) 



(48) 



These expressions are also validated by simulations. 

Since MMSE SIC with perfect decoding is equiv- 
alent to the optimum filter, we thus obtain a second 
possible equilibrium FA for the optimum filter. In 
Sec. HXl we investigate which is the FA which min- 
imizes total amount of power needed to transmit at 
equilibrium SINR. In the case of automatic ordering 
of the users, one question is whether it is best to 
order the users by increasing or decreasing total 
fading energy. The answer is the following: it is 
always best to decode the users by decreasing total 
channel energy Ei < ■ • • < Ek (see Appendix 
IXTFl) . 

An interesting feature of equilibrium FA (l47l) and 
(l48l) is that there is no limitation on the number 
of users than can be accomodated by the system, 
contrary to the previous case of (|34l) . (l38l) and (l39l) . 
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The limitation is only imposed by the increasing 
power needed for each new user decoded last, which 
grows without bound as an exponential. 

IX. Numerical Results 

In all the following, we consider that P^ax is 
chosen sufficiently high so that users can actually 
transmit at the equilibrium PA values. For the simu- 
lations, we consider the usual case of Rayleigh fad- 
ing. Although Rayleigh distribution is not bounded 
from above, simulations show that the results still 
hold. 

We consider a CDMA system with K = 32 
users and a spreading factor = 256. The noise 
variance is cr^ = 10^^°. For a number of bits 
in a CDMA packet M = 100, the goodput is 
= (i-e-^)^°° (see [24]), and jS* = 6.48. 
The capacity achieved at the Nash Equilibrium is 
C = alog2(l+/3*) = 0.39 bits/s. Unfortunately, 
the capacity itself cannot be used as a relevant 
performance measure in the definition of the utility, 
because in this case the maximal utility is obtained 
when not sending. 

We have performed simulations over 10000 re- 
alizations. Fig. \T\ shows the good fit of theoretic 
values calculated directly from (|34l), (|38]) and (|39l ) 
with those simulations. The values of the utility do 
not depend on the number of multipaths. We see 
that optimum filter requires the minimal power, and 
matched filter the maximal power to achieve the 
required goodput. 

In Fig. [2] we have plotted the average utility 
versus the number of multipaths L. Multipaths are 
supposed to be i.i.d. Rayleigh distributed with vari- 
ance 1/L, in order for the channels to have the 
same energy. Two cases are considered: the utility 
obtained in the Nash equilibrium, according to the 
VA given by dlB and ([361), and the utility in the 
case where all nodes transmit at the same power. 
For comparison purposes, the sum of the uniform 
powers is equal to the sum of the powers used 
in the Nash equilibrium. In addition, simulations 
(not reproduced here) show that this value gives the 
higher average utility for a uniform FA. The utility 
does not vary with L in the Nash equilibrium: the 
Central Limit Theorem applies to the utility, which 
is a constant times the random variable in the 
Nash equilibrium. The utility with uniform powers 
is always inferior to the utility in the Nash equilib- 
rium. However, as L increases, the gap decreases. 



as the variance of decreases, and the equilibrium 
FA becomes uniform. 

In Fig. [3] we have plotted the average of the 
inverse power of the users in the Nash equilibrium 
for each of the investigated schemes. We plot the 
average inverse power because of the direct relation 
to the utility for the users. The higher this average, 
the higher the utility for the user. The SIC filters are 
always more efficient than their linear counterparts. 
However, for a load a < 0.12 and optimum filteiB 
it is better to use the first variation of FA (|39l ) 
than use MMSE SIC (|48]). This relation is reversed 
when a > 0.12. In addition to the theoretical 
curves, Monte-Carlo simulations were performed 
both with random ordering (circles) and ordering by 
decreasing total channel energy (crosses), for L = 8 
multipaths. Simulations show that the optimal order- 
ing improves the power efficiency of the successive 
interference cancellation filters. 

In Fig. m we investigate the amelioration pro- 
vided by optimal ordering as a function of the 
number of multipaths. The simulations are done for 
K = 128 users, in order to be in the "interesting" 
zone a > 0.12. As expected, as the number of paths 
increases, the total channel energy is more and more 
the same for each channel and the gain provided 
by ordering the users decreases. However, when the 
number of users is very large and they benefit from 
automatic ordering, we see that the utility with the 
MMSE SIC equilibrium FA is the maximal utility 
that can be obtained in the non-cooperative setting. 

X. Conclusion 

Using tools of random matrices, we have de- 
rived the equilibrium power allocation in a game- 
theoretic framework applied to asymptotic CDMA 
with cyclic prefix, under frequency-selective fad- 
ing. Three receivers are considered: matched filter, 
MMSE and optimum filter (given by Shannon's 
capacity). In addition, distributed ordering mecha- 
nisms are introduced and the successive interference 
cancellation variants of the linear filters are studied. 
For each user, this power allocation depends only 
on the total energy of the channel of the user 
under consideration. For a frequency-flat channel, 
the power allocation among users is dis-uniform, 
whereas when the number of multipaths increases, 

"^The value of a is obtained as solution of the equation 
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Fig. 1. Comparison of theoretic values and simulations for utilities 
in the Nash equilibrium. 




4-1' 
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Number of Multipaths L 

Fig. 2. Simulation of utilities in the Nash equilibrium and constant 
power allocations versus L. 



Fig. 3. Average inverse power used by the different filters. 




5 10 15 20 25 30 35 



Number of Multipaths L 

Fig. 4. Simulation of utilities in the Nash equilibrium with SIC 
filter with and without optimal ordering, versus L. 



the power allocation tends more and more to a rewritten 
uniform one. 



/^(^) = / '^^'^!f,,, • (49) 



XL Appendix ° ^' Jo i+p{y) 

A. Proof of Prop. E] ^ 

2 OPT MMSE From ([B, Jq p{f,x)u{f,—a'^)df satisfies the 

^° M^M^F^^^ or^oo,^ =0,C =0 ^^j^g implicit equation dH as (3{x) and thus 
and = P{x) = 0. Thus we only have to 
prove that the derivatives of either side of (|20|) are 

equal. u{f, -a^) = . (50) 

Using p(/,x) = Pix)\h{f,x)\\ m can be /o igfef + 
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Using (|49l) and (|50l ). we can rewrite 



u{f, -a')df- 



1 



r 

Jo 



0"^ 



Jo l+/3(x)"-^ 



^ \Jo l+/3{y) +^ 



df 



-1 



P{f, x)df 

,°4Z^ + ^2 



° Jo l+/3(y) + 



+ 



dx. 



Thus from (1181 

^^OPT 



log2(e) 



-dx. (51) 



Differentiating (fT5l) with respect to a^, we obtain 

^^MMSE /-a ^ 9/3 , ,,,, 

^ = log2(e) / T—UTZK — y^)'^^- (52) 



1 



Let Tx{x) 
obtain 



l + p{x)da^ 

From dlB and ([521), we 



= - log2(e) (^P{x) + a^^(a;)^ 7r(a;)cia;. 

(53) 

From (fT3l) . we have 

° /-^ a'p{f,x)df dn 

io cr^ + jQcr^p{f,y)n{y)dyda'^ 



x)dx 



X 



' i:p{f,x)^{x)d: 

^ + lo pif^y)'^iy)dy 
1 9 /-i 



5. Influence of Other Players ' Strategies 

We want to prove that asymptotically, in the game 
{5^', P, (Mfc)fcg5A }, the strategy of a single player 
does not have any influence on the payoff of the 
other players. In other words, for all k ^ i E S^, 
for all p = (Pi, . . . , Pk) G P^, for aU G P, 



\ukip) -UkiPi,P(-i)] 

Remember that Uk 



0, as N ^ oo. 



7(A) 
Pk 



and 7 is at least 



C^. Let {Pi,. . ., Pk) be the SINRs associated with 
the power allocation p and (/JJ, . . . , /3^) the SINRs 
associated with the power allocation (P/,p(_j)). 
Then a simple Taylor expansion of 7 in P'^. gives 

im=im+{Pk-f3k)^{Pk)+o{f3k-f3k). (54) 



(55) 



According to (l54l) . it is sufficient to show that 



0, as ^ 00. 



a) Matched Filter: For the matched filter, 
the inequality is obtained directly from (fTTI) . 
The denominator of (fTTI) is always greater than 

^V"^ \h " 



I'nk 



< 



. Hence, 



|2 1, |2 



PkCT^ 



< 



P h'^ 

max "■max 



b) MMSE Filter: For the MMSE filter, the 
inequality is obtained from (fT2l) . Lemma 1 from [33] 
and Lemma 2.1 from [41], which we both reproduce 
below for convenience. 

Lemma 1: [33] Let C be a x complex 
matrix with uniformely bounded spectral radius for 
aU A^: sup^(|C|) < 00. Let w = -^[wi, . . ■,WNf 
where {wi}i=i,„N are i.i.d. complex random vari- 
ables with zero mean, unit variance and finite eighth 
moment. Then: 



log2(e) (90-2 
Observing that 



log2 I 1 



pif,y)T^iy)dy df. 



E 



w^Cw trC 

A^ 



< 



P{x) + a^^^(x) ) 7T{x)+a'^ P{x)^^{x)dx 



dn 



d 



da JO 

we obtain (|20|) from Prop. [3l 



a'^P{x)n{x)dx 



where C is a constant that does not depend on A^ 
or C. 

Lemma 2: [41] Let cr^ > q, A and B A^ x A^ 
with B Hermitian nonnegative definite, and q G C^. 
Then 

tr (((B + a^I)-! - (B + qq^ + a^I)-') A) < ^ 
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In Lemma [21 ||A|| is the spectral norm of A, i.e., 
the square root of the largest singular value of A. 
From (7), we can write 

f3k = PfcW.^Hf (G(_fc)Gf:,) + a2liv)''H,Wfe, 
P', = PfcWfc^Hf (G(_fc)'Gf^,/ + aH^y' HfcWfc 
where 



H 



G(_fe)'G(_fc) 



- G(_fc)G(^fc) 

+ (i^'-P,)(h,0w,)(h,0w, 



A corollary of Lemma [T] is that for either matrix 

C = (G(^k)Gf_i^^ + a^lN) Hfe or matrix 



C = 
[33] 



G 



'r'H ' 



Hi,, we obtain 



Wfc^Cwfc - -trC 



0, as ^ oo. 



Matrix B = G(_fc)G^^^ is Hermitian nonnega- 



k has 



tive definite, as for all w G C^, w^G 

II 1 1 2 

||G(_fc)w|| > 0. Diagonal matrix A = H^H 
spectral norm ||HfcHf || < /i^^^^. Using Lemmas [H 
and [21 as oo, we obtain 



0, as ^ oo. 



c) Optimum and Successive Interference Can- 
cellation Filters: The analog of the SINR derived 
for the optimum filter stems from the MMSE filter 
with SIC. The SINR for SIC filters have similar 
expressions with less interfering users appearing in 
the denominator. Hence the result is immediate. 



by Hence, the term P{y) \h{y)f dy in ^ 
does not depend on the actual realizations of the 
channels. Replacing P* by /3+ in ([271 ). we obtain 



that/;P(y)|%)|^rfy 



^ 



001). Replacing (3* by /3+ 
power allocation ((29l) . 



in ([281 



which gives us 
, we obtain the 



D. Expectation of the random variable ([331) 

For each user j, there are L > 1 paths with 
respective attenuations hi (-^) , £ = 1, . . . , L, which 
are i.i.d. complex random variables with mean zero 
and even distributions of the real and imaginary 
parts. The Fourier transform of those attenuations 
is = h , ^) = Eti hi ii) e-^2-^(^-i). The 
total energy of the paths is Ej = J2e=i \ he {jj) \ ■ 

We want to show that the expectation of the 



random variable X]f=i ^"^e^ equal to 1 . By 
expanding the expression of hnj, this is equivalent to 
showing that the expectation of the random variable 



hi (i) he' 



E, 



is equal to 0. Denoting by p{-) the distribution of 
hi = hi (-^), this expectation is equal to the L- 
dimensional integral of 



hihi' 



\hif + \hi'f + ^j^^ii, \hkf 



p{hi)p{hi>) Yl P{hk) 



which is an odd function of hi. Its integral is 
therefore 0, which proves the desired result. 



C. Proof of Prop. \5\ 

Given C*, we can use ([201 ) to obtain a Nash 
equilibrium power allocation in the following way. 
We rewrite ([201 ) assuming that the target SINR for 
the MMSE filter is 




= a\og,_(l + P'). (56) 

In the left-hand side of ([56[) , P{y) is given by a 
MMSE power allocation similar to the one given 



E. Proof of ^ and ([46l) 

Denote = Px-k \ hK-k\- From (|42l) . with flat 
fading, the sequence {m^jfeg^K satisfies mo = /9*cr^ 
and ruk+i = /3*cr^ + ^Z]j=o'^i- Using the fact 
that Ylt=j (j) = ( j+i) ' it immediate to prove by 
recurrence that 

Hence formula (|45l) . The demonstration is exactly 
similar for ([461 ) from the recursion tuq = (3*a'^ and 
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F. Optimal Ordering of Users 

We detemiine the ordering that makes use of the 
least total power for equilibrium PA (|45] ) (the case 
is similar for (|46l), (gT]) and (gS])). Let the ordering 

I 1 2 I 1 2 

of the users be such as \hi\ < ■ ■ ■ < . Let 
TT be any permutation of {1, . . . ,K}. Let aij = 

Then showing that the optimal ordering is such 

I 1 2 1 2 

as I /ii I < ■ ■ ■ < I /ix I is equivalent to showing that 
for any vr 

^ 1 

77— 2 '^fc^(fc) > 0- (57) 

fe=l 

Consider first a cyclic permutation. By the def- 
inition of ttij, the sum of the akTr{k) is equal to 
zero: J2k=i ^kirik) = 0. The first coefficient ai^(i) 
is positive. It is affected coefficient which is 
the greatest coefficient in the sum in (|57|) . Hence 
the sum in (|57] ) is positive in this case. 

Permutation vr can be decomposed as a product of 
disjoint permutation cycles. Each cycle determines a 
subset of indexes k, these subsets form a partition of 
{1, . . . , i^'}. With a similar reasoning as precedently, 
replacing index 1 with the smallest index in the 
cycle, the sum over the indexes k pertaining to a 
cycle of jT^akTrik) is positive. Hence the global sum 
of (1571) is also positive. 

It can be proven in a similar way that the same 
ordering maximizes the sum of inverse powers of 
the users. 
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